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A possible feature of the inert doublet model (IDM) is to provide a dark matter candidate together 
with an alteration of both direct and indirect collider constraints that allow for a heavy Higgs boson. 
We study the IDM in light of recent results from Higgs searches at the Large Hadron Collider (LHC) 
in combination with dark matter direct detection limits from the XENON experiment. We ask under 
what conditions the IDM can still accommodate a heavy Higgs boson. We find that IDM scenarios 
with a Higgs boson in the mass range 160-600 GeV are ruled out only when all experimental 
constraints are combined. For models explaining only a fraction of the dark matter the limits are 
weakened, and IDMs with a heavy Higgs are allowed. We discuss the prospects for future detection 
of such IDM scenarios in the four-lepton plus missing energy channel at the LHC. This signal can 
show up in the first year of running at </s = 14 TeV, and we present detector-level studies for a few 
benchmark models. 

PACS numbers: 14.80.Cp, 95.35.+d 



I. INTRODUCTION 



The era for studying particle physics with the LHC at 
CERN is ongoing. Since 2010, the experiments have been 
collecting data from proton-proton collisions at a center- 
of-mass energy of 7 TeV. This has already enabled the 
exploration of new regimes of the current standard model 
(SM) as well as physics beyond the SM. One of the aims 
is to establish or exclude the presence of a SM Higgs 
boson 1 . The latest public Higgs search results were pre- 
sented bythe ATLAS and CMS collaborations in March 
2012 [HHf. These analyses exclude a SM Higgs boson 



'Electronic address: mgustafs@ulb.ac. be 
t Electronic address: sara.rydbcck@desy.dc 
tElectronic address: llopezho@mpi-hd.mpg.de 
§ Former academic affiliation. 

The term Higgs boson will throughout the text be used for the 
physical scalar particle emerging from the electroweak symmetry 
breaking in the SM by the Brout-Englert-Higgs mechanism [F. En- 
glert and R. Brout, Phys. Rev. Lett. 13, 321 (1964); P.W. Higgs, 
Phys. Rev. 145, 1156 (1966)]. 



in the range 127 — 600 GeV to the 95% confidence level 
(CL). It is however important to keep in mind that new 
particles can both contribute to the Higgs decay width 
and alter its production cross section. The exclusion lim- 
its on this range of Higgs boson masses might thus not be 
valid for a Higgs boson that is SM-like in many respects, 
but which also couples to states beyond the SM. This is 
of particular relevance to the present paper. Let us point 
out that while both the ATLAS and CMS experiments 
have started to see potential evidence for a particle signal 
at ~ 125 GeV, the significance is not yet enough to claim 
discovery and establish this to be caused by the Higgs 
particle itself. Moreover, there have been other, perhaps 
interesting, excesses in the Higgs searches; e.g. at the 2a 
level for a ~ 320 GeV particle mass in CMS [H, H[ , which 
was then not confirmed by the latest preliminary results 
from the ATLAS experiment 

Because of the nature of hadron colliders, the LHC has 
obvious advantages in probing beyond the SM scenarios 
that incorporate strong quantum chromodynamic (QCD) 
interactions, such as minimal low-energy supersymmetry 
models. So far the LHC searches have found no evidence 
for strongly interacting beyond the SM particles @ . No- 
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toriously, scenarios without direct SM QCD interactions 
are expected to give lower signals - although many ex- 
ceptions, such as resonances (e.g. Q) or composite state 
effects (e.g. @,@), may appear. From an empirical point 
of view, there is a priori no need for new QCD inter- 
acting sectors. Indeed, two of the major questions in 
particle physics and cosmology - the fine-tuning problem 
in the SM Higgs sector (commonly known as the "LEP 
paradox" or the hierarchy problem [l(|), and the dark 
matter (DM) problem with a thermally produced weakly 
interacting massive particle (WIMP) as one of the long- 
standing candidate solutions jLlT - [l5l ] - are not directly 
connected to QCD properties. 

Given the latter point of view, we study the inert dou- 
blet model (IDM) ; a minimal extension of the SM which 
contains one additional electroweak scalar doublet and 
has the potential both to alleviate the mentioned fine- 
tuning in the SM and to provide a DM candidate. The 
IDM appeared already in the 1970s i n (lo| . but received 
more attention after Barbieri et.al. (see also (l8j ) 

showed that the model could provide both a DM can- 
didate by an imposed Zi symmetry and allow for SM- 
like Higgs masses up to 600 GeV without contradicting 
electroweak precision data. These authors pointed out 
how raising the Higgs mass could alleviate the problem 
posed by the LEP paradox [ljj and thus eliminate the 
fine-tuning in the SM up to an energy scale of a few 
TeV (see however [lj|). Regarding its DM candidate 
[20h22| . many signatures have been studied and range 
from potentially striking gamma-ray lines plj . cosmic- 
ray and neutrino fluxes [23|, [24| to direct detection sig- 
nals [H E! El- The lack 01 conclusive beyond the SM 
signals in these channels and in the data from collider 
experiments [13, HH, H3, HH has so far only partially con- 
strained the IDM parameter space. 

We devote the first part of this paper to the question 
of whether a large Higgs boson mass (> 160 GeV) can 
still be accommodated within the IDM. Indeed, fairly 
large inert particles-Higgs couplings are needed in order 
to elude the SM Higgs searches at the LHC. The same 
couplings are, however, severely constrained by DM di- 
rect detection searches at XENON-100. 

The need for large scalar couplings leads us to the sec- 
ond part of the paper, where we study a new potential 
discovery channel for the IDM in the form of multilcpton 
events via heavy Higgs production at the LHC. 

Even if the picture of the Higgs being SM-likc will be- 
come clearer as the LHC continues to run during 2012, 
the possibility and the nature of a modified Higgs sector 
might remain an open question. After 2012, the LHC 
will have a long shutdown in preparation for start-up in 
late 2014 at the design center-of-mass energy of 14 TeV. 
Once the design luminosity is reached, the experiment 
will accumulate ~ 100 fb _1 per year, allowing probing of 
the Higgs sector and beyond SM physics in more detail. 

The prospects for detecting IDM signatures in the up- 
coming LHC data at 14 TeV has already been partly ex- 
plored. In [13, H3] the authors studied how inert particles 



affect SM Higgs searches, by the opening of additional 
decay channels, as well as the discovery potential in the 
dilepton and missing energy channel. A more compre- 
hensive study of this dilepton channel was done in [29J, 
followed by a trilepton study [3(| ■ None of these studies 
explore the possibility to detect the inert doublet model 
in the almost background-free multilcpton (> 4 leptons) 
plus f/T channel. Here we argue that it is natural to 
study the tetralepton channel in addition to the dilep- 
ton and trilepton channels. This has actually been done 
for many other popular models, e.g. in supersymmetry 
[3T| - [35l | and extra dimension [H, H3] models. 

The inert doublet contains four new particle states. 
The more massive states may be pair produced in pro- 
ton collisions and subsequently cascade decay (in one or 
two steps) down to the lightest inert particle state, which 
remains stable due to the conserved Z2 parity. In each 
decay step, an electroweak gauge boson is produced and 
can decay into one or two charged leptons. If the lightest 
stable inert particle is electrically neutral, it will con- 
tribute to the missing transverse energy (J£t)j and up to 
six charged leptons can be directly produced from the 
W ± and Z boson that participated in the cascade decay. 
We show that the (> 4/ + ~$t) channel is an interest- 
ing test of the IDM and can provide an early discovery 
channel of the IDM when the LHC runs at 14 TeV. 

In Sec.HIland lllll we set up the IDM framework and the 
theoretical, experimental and observational constraints 
that will be imposed on the model. In Sec. IIVI we an- 
swer our first question, namely under what conditions 
a heavy SM-likc Higgs can survive the recent and com- 
plementary constraints from the LHC and XENON. In 
Sec. [V] we turn to our second aim, to discuss the multi- 
lcpton signal at the LHC in such scenarios. We perform 
detailed event simulations for a set of IDM benchmark 
models and the SM background and describe our analy- 
sis tools in Sec. I VI I Our results and discovery prospects 
for IDM in the tetralepton+^?T channel are presented in 
Sec. I VI El and in Sec. IVIII we summarize and conclude. 



II. THE INERT DOUBLET MODEL 

The IDM consists of the SM, including the standard 
Higgs doublet Hi, and an additional Lorentz scalar in 
the form of an SU(2)l doublet if 2. An extra unbroken 
Z2 symmetry is introduced, under which H2 is taken to 
be odd (H2 — > —if 2) while ifi and all other SM fields are 
even. This Z2 symmetry protects against the introduc- 
tion of new flavor changing neutral currents and guaran- 
tees the absence of direct Yukawa couplings between the 
inert states and the SM fermions (hence the name inert 
doublet model). The symmetry also renders the light- 
est particle state of if2 stable. If neutral, the latter can 
provide a good DM candidate. The new kinetic gauge 
term takes the usual form, D^H^D^H^, and the most 
general renormalizable CP conserving potential for the 
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IDM scalar sector is 

V = ^\H 1 \ 2 + f x 2 2 \H 2 \ 2 + X 1 \H 1 \ 4 + X 2 \H 2 \ 4 

+ X^mflHil 2 + \i\HlH2\ 2 + \ 5 Re[(HlH 2 )% (1) 

where /i 2 and A^ are real parameters. 

Four new physical particle states are obtained in this 
model: two charged states, H ± , and two neutral states, 
H° and A . After standard electroweak symmetry break- 
ing, the masses of the scalar particles (including the SM- 
likc Higgs mass m^) are given by: 

m 2 H0 = [i 2 + (A 3 + A4 + A 5 )u 2 = [il + \ H av 2 , 



m A° — M2 + (A3 + A4 - A 5 )v 2 = /ij + Xa«v 2 



m 2 H± 



M 2 + , 
-2n 2 = 4A lW 2 , 



(2) 



where v sa 177 GeV is the vacuum expectation value 
of the Higgs field H\. In the following, we choose H° 
to be the lightest inert particle, and hence the potential 
DM candidate. Notice that the roles of ^4° and H° are 
equivalent in the IDM and our conclusions would remain 
unchanged if we had chosen A to be the DM candidate. 
A convenient set of parameters to describe the full scalar 
sector are the four scalar masses {m^o, to^o, m H ± , m^}, 
the self-coupling A2 and A#o = A3 + A4 + A5 . 



III. CONSTRAINTS ON IDM 

There are several theoretical, experimental and obser- 
vational constraints on the model that have to be con- 
sidered. For all the models in this study we consistently 
impose: 

• the requirements for vacuum stability (38l l39j , 

• that calculations should be within the perturbative 
regime (with A, < 4tt) [13, HI] 2 , 

• unitarity constraints (the absolute value of the eigen- 
values of the S matrix are required to be < 1/2 for 
scalar-to-scalar scatterings, including the longitudi- 
nal parts of the gauge bosons) |40l445| 3 . 

• consistency with electroweak precision tests 
(EWPT) (99% CL) H3, 

• consistency with particle collider data from LEP 

(-95% cl) [Tsui mm, 



2 The constraint in Eq. 17 of reference llTjl . that poses a suffi- 
cient condition not to affect their naturalness arguments for the 
IDM, is not included. Applying it does not change our conclu- 
sions, although it would reject the models in our scans which have 
m H° ~ 120 GeV and correct relic density Q H o ~ SJcdm- 

3 See also [46ll . where the authors studied the constraints from uni- 
tarity on the IDM. 



• a relic abundance of H° in agreement with the 
WMAP measured f> C DM^ 2 = 0.1109 ± 0.017 (3a) 

H, 

• consistency with direct DM searches by XENON 

(90% CL) |3, mi, 

• consistency with indirect DM searches. We in- 
clude the 95% CL gamma-ray constraints by Fermi- 
LAT (assuming Navarro-Frenk- White profiles) (52r - 
|5"H - No other indirect detection probes are consid- 
ered here as these either give significantly weaker 
limits or are associated with too large astrophysical 
uncertainties. 

IDMs with large Higgs masses can potentially alleviate 
the fine-tunin g pr esent in the SM and thus address the 



LEP paradox [17J- While we choose not to impose any 



explicit naturalness constraints here, we will extensively 
comment on this in Appendix [AJ 

For a review of many of the constraints on IDM we 
refer to [39[ . We have implemented the above list of con- 
straints, as described in the given references, into our own 
computer code. We stress the importance of combining 
all these bounds since, as we will see, their complemen- 
tarity becomes a powerful tool in constraining the IDM. 

We will present results of random scans over the whole 
viable IDM parameter space that is of interest for our 
study (from a few GeV to hundreds of GeV). More pre- 
cisely, the free parameters were taken to be the three 
masses of the inert scalars, the Higgs mass and the cou- 
pling Al. We scanned over the ranges: 

5 GeV < m H o < 170 GeV, 

m H o < m A o < 800 GeV, 

max(m ff o, 70 GeV) < m H ± < 800 GeV, 

100 GeV < m h < 900 GeV 

10~ 5 < |A L | <4tt. 

Once these parameters were chosen randomly, the value 
of A2 was fixed to its minimal value satisfying the con- 
straints from vacuum stability. The resulting IDMs were 
confronted with the constraints listed in this section and 
only models passing the full set of constraints were con- 
sidered as viable. 

A random scan is always incomplete in covering all pos- 
sible models. By a combination of random scans, simple 
Markov chain Monte Carlo searches (following [531) an< ^ 
physical insight into where models could be expected to 
be found, we believe that we have been able to cover all 
relevant parts of the parameter space for our results with 
more than 100 000 models. For example, earlier studies 
[56l . |60| have already shown that expanding the scan to 
larger H° masses is not relevant if H° should constitute a 
WIMP DM candidate. This is at least true for Ho masses 
below 500 GeV, and higher masses are not relevant for 
the current LHC searches. It is worth noting that this 
part of the IDM gives well isolated regions in all our pre- 
sented quantities. In practice, no viable IDMs were found 
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with m h > 700 GeV, m H o > 150 GeV, m A o < 50 GeV 
or |A H o| > 7. 

The DM relic density calculations have been performed 
by DarkSUSY [HtJ interfaced with FormCalc 53- This 



Dark matter direct detection limits 



code was originally developed in [2lj . but has now been 
updated to also include three-body final states (as in 
|55j). Also, an upgrade of micrOMEGAs [59] including 
annihilation into three-body final states [55j has been 
used for the scans. 



IV. IDM IN LIGHT OF XENON AND THE LHC 
HIGGS SEARCH 

Dark matter direct detection and the LHC's SM Higgs 
searches are known to be complementary in constraining 
Higgs portal DM models [6ll - l69| . Direct detection ex- 
periments pose upper limits on the DM coupling to the 
Higgs. This in turn restricts the Higgs decay rate into 
the invisible DM states, which makes it more difficult for 
such models to escape the bounds coming from LHC's 
Higgs particle searches. 

The constraints on singlet scalar DM from combining 
XENON-100 and the LHC SM Higgs searches were e.g. 
studied in [(32| for a wide range of Higgs masses. Let us 
emphasize that the latter analysis did not assume any ex- 
plicit mechanism for evading EWPT constraints, which 
would otherwise constrain the SM Higgs mass to be be- 
low roughly 160 GeV. By contrast, the IDM provides such 
a mechanism and can easily accommodate Higgs masses 
up to at least 600 GeV while still being in agreement 
with EWPT. Another difference to the singlet scalar DM 
model is that the IDM's "dark" sector is composed of 
more than one particle state. The additional states po- 
tentially provide new contributions to the decay width of 
the SM-like Higgs boson, along with additional processes 
relevant for the determination of the DM relic density. 



A. Constraints from direct detection DM searches 

Figure Q] shows how XENON [H3, HH constrains the 
IDM models that have a relic density in agreement with 
WMAP. These constraints assume a local H° density of 
po = 0.3 GeV/cm -3 and a standard Maxwellian veloc- 
ity distribution. The spin-independent cross section for 
IDMs is calculated as in fl7l: 



J5I 

a H°-p 



_ ; (3) 
47t(to„ + m H o) 2 m\ ' 

where the form factor is taken to be / = 0.3 [20l. 1(321 Pzol] . 
and m n is the target nucleon mass. The loop induced 
contribution estimated in fl7j is also included, but it is 
very small. 

This leaves a viable mass range roughly between 45 
and 80 GeV for the DM candidate H°. The range 
can be extended up to ~ 150 GeV with a few mod- 
els marginally surviving the current XENON-100 bound 




m DM [GeV] 



FIG. 1: Direct detection signal for IDMs in agreement with 
WMAP data. The crosses represent IDM models that pass 
all our imposed constraints from Sec. IIIII (not including LHC 
constraints). The upper 95 % CL bounds from XENON-100 
[U and XENON-10 [H are shown by the solid lines. The 
points labeled A1-C4 show the scattering cross section for our 
benchmark models. As explained in Sec. lVIl the B models can 
only pass the constraints after taking systematic uncertainties 
into account, and the C models account only for a fraction of 
the DM density. 



[56j. The low mass region below < 10 GeV is excluded 
both by XENON-10 [B(J and by Fermi-LAT gamma-ray 
constraints [H, [7l|. 4 We will, however, include low H° 
masses in parts of the following discussion for illustrative 
purposes, although they are excluded once we impose all 
our constraints. 

A viable large H° mass region above ^500 GeV also 
exists [60], but is not of interest for the present study. 
Such heavy IDM states would for kincmatical reasons 
never alter the width of the Higgs boson (with a mass 
below 1 TeV) and therefore the LHC constraints apply 
exactly as in the SM. Such heavy IDM states will also 
be very difficult to probe directly at the LHC. On top 



4 Concerning a low mass WIMP, there is a debate as to what extent 
the exclusion limits from direct detection results are reliable (see 
e.g. [73n . In order to be conservative, we could therefore choose 
not to include the XENON-10 upper bounds. At the same time, we 
note that the WIMP signal constraints from the Fermi-LAT data 
on gamma-rays from e.g. dwarf galaxies [53i, l7ll] also exclude this 
low H° mass region of the IDM. We therefore take the viewpoint 
that a light H° below 10 GeV is not a viable WIMP candidate 
within the current standard scenario [391 . 



5 



of that, in order to get the correct relic density and to 
comply with EWPT, only small Higgs masses can be con- 
sidered 1601. 



B. Constraints from Higgs boson searches 

The latest results are based on analyses of ~ 5 fb _1 
of integrated luminosity. The CMS experiment set the 
strongest (preliminary) constraints on large Higgs masses 
until March 2012, excluding a SM Higgs boson over the 
mass range 127-600 GeV to 95% CL, when all search 
channels are combined (4.6 — 4.7 fb _1 of integrated lumi- 
nosity) [11 • At that time, ATLAS presented their (pre- 
liminary) limits on large Higgs masses using up to 4.9 
fb _1 [lj]. The CMS collaboration also updated their lim- 
its in some channels for 4.6 — 4.8 fb _1 [2[. We will here 
use both the experiments' current best exclusion limits 
on a Higgs signal ct/ctsm- Here (t/ctsm denotes the signal 
rate in units of the expected SM Higgs production cross 
section osm- The 95% CL upper limits, for all channels 
combined but for each experiment individually, will be 
used. In Figure [2] the excluded signal strength, as a func- 
tion of the Higgs boson mass, is shown as the blue (gray) 
region. The exclusion region represents the strongest of 
the two limits from the CMS (dotted line) and the AT- 
LAS (dashed line) experiments. 



1. Reduction of the Higgs signal in the IDM 

In the IDM, the new contributions to the SM-like Higgs 
width F/i can have a significant impact on the LHC 
Higgs searches by effectively reducing the Higgs produc- 
tion cross section into SM particles. Since H° is neutral 
and stable, the Higgs decays into H° pairs will neces- 
sarily contribute to an invisible width. However, let us 
emphasize that the Higgs can also decay into A and H ± 
pairs which would further increase the Higgs width. The 
latter processes give rise to the production of (off- or on- 
shell) Z and W bosons that can make them partly visible 
in the Higgs search channels. 

Nevertheless, the exclusion limits on the Higgs mass 
range could very well be evaded within the IDM. The 
processes h -> H°H°, h -> A A and h -> H+H~ en- 
hance the Higgs decay width by: 



for producing SM particles xx is given by 



AF 



IDM 



16Trmh 



4m\ a 



V 



1/2 



4m 2 H „ 



1/2 



2Xi 1 - 



4m 



i/2- 



,(4) 



ctidm(w -> h) Br(h -> xx)i DM 



&sm(pp -> h) Bi(h 



SM 



F 



SM 



XX) S M 



(5) 



F 



SM ' 



Ar 



IDM 



pSM 

h— >xx 



where the sum runs over <j> = A°,H ± and TgM is the 
total decay width of the Higgs in the SM. In the last 
step we use the fact that the Higgs production and decay 
rates into SM particles xx are unchanged to first order 5 
(r™ SM = r J ™ SM ). e denotes the efficiency with which 
A A and H + H~ may contribute to the current xx Higgs 
search. This efficiency may be expected to be low due to 
the fact that the final states will contain extra invisible 
H° states and therefore, in principle, have different char- 
acteristics than the pure SM xx final states. 

This means that for the whole range of Higgs masses, 
even if excluded within the SM, the LHC limits could 
potentially be evaded within the IDM. 

In the next subsection we will argue that, for the mod- 
els of interest for our study, Higgs decay into all IDM 
particles will effectively be invisible. In that case the 
reduction factor in Eq. (J5j> reduces to 



SM 



F 



SM 



Ar 



IDM 



(6) 



In general, the IDM contribution could, in principle, also 
enhance certain SM Higgs signatures, depending on the 
specific model and search channel. However, the effect 
of such a contribution would only give stronger exclusion 
limits on large Higgs masses than in the SM. Taking TZ = 
JZ cons is thus the most conservative choice when it comes 
to determining to what extent the IDM is excluded and 
therefore the one that we will adopt in the following. 



2. Higgs searches and the IDM 

The WW and ZZ search channels are the most ef- 
fective ones in the search for heavy Higgs bosons, and 
below we list their most sensitive subchannels. We quote 
the excluded SM Higgs masses, as this indicates where 
the searches could be sensitive enough to exclude Higgs 
masses in the IDM. 

• h -> ZZ^ — s- 41 with I = {electron, muon}. By 
these lepton channels alone, the CMS experiment 
excluded at 95% CL SM Higgs boson masses in 
the ranges 134-158, 180-305 and 340-465 GeV 0. 
At the same confidence level, ATLAS excluded the 
ranges 134-156, 182-233, 256-265 and 268-415 GeV 



where Xh° A°,3 are given in Eq. (JXJ) and (|2|). 

In the narrow- width approximation, the signal 
strength ct/ctsm, or equivalently the reduction factor 1Z, 



5 For an effect at the loop level, see e.g. the study in [73j of the 77 
channel. 
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Upper/lower limits on the Higgs signal in IDM for different H masses 




■ Excluded by LHC Higgs searches 
CMS (dotted), ATLAS (dashed) 



- Lower limits from XENON (IT 100% of local DM) 
for m u o = 8, 45, 150 GeV 



150 200 250 



300 350 400 
m h [GeV] 



450 500 550 600 



Upper/lower limits for m H o = 70 GeV, under different assumptions 
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[Excluded by LHC Higgs searches: 
CMS (dotted), ATLAS (dashed) 
Lower limits from XENON for m H o = 70 GeV 
(with H° 1 00% of the local DM) (a) 
■XENON limits shifted by syst. uncert. (b) 
■ If H° constitutes 1 0% of the local DM (c) 



150 200 250 300 



400 450 500 550 600 



eV] 



FIG. 2: Left: The solid lines show, for different H° masses, the lowest limits on the Higgs production rates a/asM that 
the IDM can have at the LHC and still be compatible with XENON- 100 data. These limits correspond to when none of 
the produced IDM states contribute to a signal in the LHC's Higgs boson searches, i.e. 1Z = TZ cons . The constraint from a 
thermal freeze-out calculation of H° has not been applied yet, and H° is assumed to make up the local DM density. All other 
constraints, EWPT and LEP limits in particular, are taken into account. The dark (blue) region is excluded by the current 
Higgs search results from the LHC. Thus, only the white regions above the lines remain allowed in the IDM. Right: All lines 
assume an H° mass of 70 GeV. The upper (red) line assumes that H° particles provide all the local DM, the dark (blue) 
line below is after systematic uncertainties (see text) are included to diminish the direct detection constraints, while the lower 
(green) line applies when assuming that H° contribute only 10% to the local DM density. 



Q . An important requirement in all these searches 
is that at least one same- flavor opposite-sign lepton 
pair has an invariant mass in a window around the 
Z mass. 

• h —> ZZ —> 2l2v. In this channel, the events are 
required to contain a minimum amount of missing 
transverse energy and the lepton pair is required to 
form an on-shell Z boson. The CMS collaboration 
was able to use this channel alone to exclude SM 
Hig gs m asses in the range 270-440 GeV at 95 % 
CL [74j . and the corresponding range excluded by 
ATLAS is 320-560 GeV [z|. 

• h —> ZZ —> 2l2q. This search requires the invariant 
mass of the jet pair to correspond to an on-shell 
Z boson. ATLAS excludes SM Higgs masses in 
the ranges 300-310 and 360-400 GeV at the 95% 
CL [Z6[, while CMS could not exclude SM Higgs 
production cross sections by the use of this channel 
alone fn\ . 

• h —> WW —> Ivlv. In this channel, the events are 
required to contain at least two leptons of oppo- 
site sign and missing transverse energy. Cuts on 
the transverse mass, reconstructed from the lepton 
pair together with the missing transverse energy, 
are also applied. The CMS excludes at 95% CL a 
SM Higgs mass in the range 129-270 GeV Q, and 
ATLAS the range 130-260 GeV 0. 



Higgs decays into H° would be invisible, but might 
it also be the case that decays into A and i7 ± escape 
detection in the above search channels? 

The decay channel h A A would give rise to two 
Z bosons and could be visible in the above ZZ search 
channels. It would however only give a visible contribu- 
tion if (m^o — m H o ) is large enough to produce on-shell 
Z bosons via the decay A ->• H° + Z. 6 

In the WW — > 2l2v Higgs search channel, the final 
state is required to include two opposite-sign leptons and 
missing energy. The h — > A°A° and h — > H + H~ pro- 
duction, with the subsequent decays A — > H° + Z and 
i? ± — > H° + W ± , could pass these requirements and one 
can imagine that this could contribute to a signal in this 
search channel. Let us therefore take a closer look at this 
possiblity, to see if the contribution could be significant. 
So far this channel only excludes SM Higgs masses in the 
range 130-270 GeV, and we therefore expect that it is 
only within this same mass range that Higgs bosons can 
be excluded in the IDM. This statement is motivated by 
the use of cuts on the transverse mass, that sets the SM 
Higgs mass for which the limit applies. This 'transverse 



Even in the case of on-shell Z's, the characteristics of the final states 
are altered by the presence of H°'s giving rise to In the ZZ — > 
41 channel this would lead to a smearing of the 4Z invariant mass 
spectrum, thereby evading a peak search, but could potentially 
contribute to the observation of a less constraining broad excess. 
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mass' variable corresponds to the Higgs boson mass in 
the SM and should roughly do so also in the IDM. This 
entitles the use of the same a I asm limit for the Higgs in 
the IDM as in the SM. 

In this specific mass range, Higgs decays into A A and 
H^H^ will however never contribute to the WW — > 2l2v 
Higgs search channel. This is because for < 160 GeV 
the LEP limits [28], already exclude almost all in- 
ert particles A and H ± with masses less than 80 GeV, 
which are the only masses that could have been kinemat- 
ically accessible for these Higgs decays. The exception, 
with lighter m^o #o, occurs only when the mass splitting 
m A a — m H a is very small, and the final-state fermions 
are then too soft to contribute. Moreover, in the re- 
gion 160 GeV < m h < 270 GeV it turns out that IDMs 
which account for all the DM are excluded irrespective of 
whether the A and H ± states are invisible to the Higgs 
searches or not (see Figure [2]). 

This means that for many models, in particular those 
that have a mass difference (m^o — mjjo) too small to 
produce Z bosons on shell, the IDM contributions to the 
Higgs width can be treated as invisible in the current 
LHC searches for heavy Higgs bosons. Our arguments 
for such a treatment were based on the channels impor- 
tant for the searches in the high mh region, while for low 
Higgs masses, other channels could be more important. 
Nevertheless, we will apply the same assumption to all 
our models as this will not alter our discussion. 
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FIG. 3: All the points correspond to IDMs giving rise to a DM 
relic density in agreement with WMAP. Small (grey) points 
represent models that neither pass the LHC Higgs searches 
nor XENON- 100 constraints, diamonds (in blue) represent 
models that pass LHC and large (green) points represent mod- 
els that pass XENON-100 constraint individually. Squares (in 
red) represent models that pass the constraints from both the 
LHC and XENON-100 experiments. All the other constraints 
of Sec. IIIII have also been taken into account. 



C. Constraints on IDM from LHC and 
XENON-100 combined 

Figure [5] shows the LHC Higgs exclusion limit together 
with IDMs that have the largest invisible Higgs width 
possible and still pass XENON-100 direct detection con- 
straints. As we can conclude from the above discussion, 
all the inert states resulting from Higgs decay can be re- 
garded as effectively invisible when the mass difference 
(m^o — m H o) is less than niz, i.e. 1Z = lZ cons . In Fig- 
ure [U we present lines for when we take 1Z = JZ cous for 
some representative to#o masses. 

Once m#o and m>h are fixed, Eq. ([3]) and the XENON- 
100 exclusion limit on cr SI determine the largest available 
value of Ajjo, and consequently Ar^^o^o . The largest 
values of A^o and A3, driving the two other contribu- 
tions to AriDM in Eq. [4] can be found numerically under 
the imposition of all the other IDM constraints listed in 
Sec. IIIII The only exception is that we do not yet impose 
that H° accounts for the total WMAP DM relic abun- 
dance. Instead, wc immediately assume that the local H° 
density provides the observed DM density that is relevant 
for the constraints on DM direct and indirect detection. 
This is in order to keep the discussion more general at 
this stage, and not include constraints from the freeze-out 
process occurring in the early Universe. We notice that 
the LEP and EWPT bounds give the most crucial lim- 
its to constrain A^o and A3 after the XENON bound has 



been imposed. Together with the XENON and LHC con- 
straints, they are efficient in excluding IDMs with heavy 
Higgs masses. 

We see that even without including the relic density 
calculations, the XENON and LHC Higgs searches, if 
taken at face value, exclude most of the IDM scenario 
with large Higgs masses. Only two exceptions appear - 
see the left plot in Figure [5] 

First, we have the low mass WIMP, with e.g. mjjo = 8 
GeV, which could give rise to large Higgs decay branching 
ratio into H°. As discussed in Sec. IIV Al this case is al- 
ready excluded by XENON- 10 and Fermi-LAT data and 
is presented for illustration only. The second exception 
arises in the large mass region for m H o ~ 80 — 150 GeV, 
which might still be viable for the largest Higgs masses. 
However, if we take into account also the constraint from 
having the DM candidate, H°, as a thermal relic, this re- 
gion is no longer allowed. This is clearly seen in Figure [3] 
where the relic density calculation has been included. We 
are thus able to exclude the so-called 'new viable region' 
of IDM found in [56| even before direct detection exper- 
iments have fully probed this regime of the IDM. There- 
fore none of these exceptions provides good models. 

Also the possibility to have models with Higgs masses 
above 600 GeV still remains. The LHC has only pre- 
sented bounds on er/ersM for Higgs masses below 600 
GeV, and we can therefore not use this quantity di- 
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rectly to exclude models with such large Higgs masses. 
The EWPT and unitarity constraints, however, limit 
the Higgs mass to be below ~ 700 GeV (also the triv- 
iality /perturbativity bound would disfavor larger Higgs 
masses [13,[8l|). As can be seen from Figure [3l when the 
thermal relic density calculation has been included, the 
DM mass range m^a ^45 — 80 GeV with a very heavy 
Higgs in the range 600-700 GeV is still an allowed region. 

In Figure [3] we present the result of a random scan in 
the mom € [15 — 170] GeV parameter space of the IDM 
giving rise to an H° relic abundance in agreement with 
WMAP HI] at the 3cr level. All the constraints from 
Sec. IIIII are now included. The plot illustrates in the 
mh — m H o plane the IDMs that pass the constraints set 
by XENON-100, the LHC Higgs searches and WMAP. 
We see that many models pass either the direct detec- 
tion or the LHC Higgs bounds individually. In the heavy 
Higgs region there are no surviving models, except for 
the region m/j > 600 GeV and m#o ~ 45 — 80 GeV (see 
also the plots in Figure [2]) . We thus conclude that in or- 
der to have an IDM that makes up all the DM and has a 
SM-likc Higgs boson in the 160-600 GeV mass range, at 
least one of our imposed constraints has to be relaxed. 



D. Accommodating a heavy Higgs boson 
and DM in the IDM 

One of the original motivations for studying the IDM 
was that it could alleviate the LEP paradox in the SM 
by allowing for a heavier Higgs particle while staying in 
agreement with EWPT. We have shown above that con- 
straints from direct detection in combination with the 
SM Higgs search essentially rule out large Higgs masses 
up to - 600 GeV in the IDM. 

In this section, we investigate the assumptions that 
could be relaxed in order to allow for a large range of 
high Higgs masses (mu > 160 GeV) within the IDM. 
In particular, we will allow for larger values of Xh° by 
suppressing the bound that derives from direct detec- 
tion searches. In that way, models with larger invisible 
Higgs branching ratios will become available, which con- 
sequently give lower signal strengths in the LHC Higgs 
searches. This is illustrated in the right panel of Figure [5] 

The bound on Ajjo can be suppressed in two ways: 

1. by assuming that the DM from IDM docs not ac- 
count for the entire DM abundance: the green line 
in the right panel of Figure [2] assumes that H° con- 
stitutes only 10% of the local DM density po- This 
suppresses the constraint on cr^ _ p by the same 
factor. 

2. by considering systematic uncertainties in direct 
detection: the dark blue lines in the right panel 
of Figure [5] take into account a smaller form fac- 
tor / = 0.26 [H, IH,_a smaller local DM density 
po = 0.2 GcV/cm -3 [83||, and, in addition, include 



a minor effect of 10% weakening of the XENON- 
100 cross section limits due to uncertainties in the 
local WIMP velocity distribution [Hj]. Concerning 
the local DM density, there have been recent im- 
proved measurements constraining it to the range 
po = 0.3 ± 0.1 GeV/cm" 3 (see also (H). 

These reconsiderations weaken the constraints on A#o , 
and IDMs with irih > 500 GeV could be allowed. Also 
to/i around 320 GeV could be allowed if only the LHC 
constraints from CMS are considered. However, the pre- 
liminary analysis recently presented by ATLAS [l[ does 
not show any excess around rrih = 320 GeV as CMS 
does, but instead puts very strong constraints in the 300 
to 450 GeV mass range. There are also uncertainties re- 
lated to the absolute calibration of cross section limits 
at the LHC on d/osM- We choose here not to take into 
account such potential additional uncertainties. 

If H° particles constitute only a fraction of the DM 
density they would more easily pass direct detection con- 
straints (now rescaled by ^Icdm/^h ) while having a 
larger A#o coupling, and then be able to evade the LHC 
Higgs limits. It then remains to be shown if such models 
exist that have such a low relic density while not ex- 
ceeding the other constraints in Sec. IIIII The possible 
mechanisms for this in the IDM are as follows 7 : 

• Annihilation via h at the resonance (mno ~ m^/2): 
In the case of a heavy Higgs boson, the resonance 
could only occur when m H o > 80 GeV and an- 
nihilations into gauge bosons already provide an 
efficient annihilation mechanism. 

• Coannihilations {m H o ~ m^o or m H a ~ m H ± ): 
This is relevant for small mass differences when 
m A°,H ± / rn H° ;$ 1-1. For large Higgs masses, the 
EWPT also requires that (rriH± — m^o) x (m^o — 
m H o ) is positive fl7j . This means that m H ± > m^a 
and that the mass difference between the two neu- 
tral inert scalars has to be small. For the tetralcp- 
ton search channel that wc will investigate in the 
next section, this has the implication that the lep- 
tons from the decay A° — > H Q are too soft to be 
detected at the LHC. 

• Annihilation to WW, ZZ and tt {m H o > m^): 
Strong annihilation channels into gauge bosons be- 
come kinematically available already for m#o just 
below mw, nnz or m t . 



7 Models with annihilation dominantly into fermions have (<jv) oc 
A^o , and are already in the region excluded by direct detection 
searches. This can be seen in Figure [T] where m H o < 40 GeV cor- 
responds to models having annihilations into fermions only. In that 
framework, increasing \ H o would not alter the bounds from direct 
detection searches. Indeed, these bounds derive from the quan- 
tity a SI xfl H o oc \ 2 H0 /(av) which is unchanged under a rescaling 
of A^o . 
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Although all three of the above mechanisms could be 
viable, we will in the next section only consider models 
where the relic density is suppressed by the last type of 
mechanism. This is because we want to investigate the 
best prospects for detecting the IDM in the tetralepton 
channel at the LHC, and the simplest scenario to consider 
is then when the WW annihilation channel regulates the 
DM abundance. 

We will also consider benchmark models that give a 
relic density in agreement with 100 % of the observed 
DM. However, for these models systematic uncertainties 
for the direct detection searches have to be included, as 
described above, to make them pass all constraints. 

V. THE MULTILEPTON SIGNAL 

The inert scalars can only be produced in pairs, since 
each inert particle has negative Z2 parity contrary to the 
SM particles. At tree-level, the relevant hard processes 
producing final states with four leptons or more, are via 
the gauge bosons and the Higgs boson: 

qq' -> W ± -> A°iJ ± (7) 
qq' -> Z/j/h -> H + H~ . (8) 

The tree-level contribution to qq' — >• h — > A A is negligi- 
ble but at loop level, gluon fusion into Higgs is important 
for A A and H + H~ production. 

After the inert particles are produced, they will cas- 
cade decay through the processes: 

H±^{f Q ™t , and A°^H°Z (9) 

or 

H ± H o w ± t and A Q ^ I H±WT ^ (iQ) 

depending on whether or ^4° is the most massive 
inert state. The gauge bosons will, with their respective 
branching ratios, decay into fermions, /, according to 

W ± -> / ± i/ and f ± f ± . (11) 

Figure |4] illustrates these production and decay chains. 
Our focus will be on the production of four or more lep- 
tons, I (which in this context refers only to electrons and 
muons), where the SM background is expected to be very 
low. 

The cross sections and decays widths will be calcu- 
lated using MadGraph/MadEvent [86j |. and their stream- 
lined interface with Pythia [I?} and PGS [88[ to simu- 
late hadronization and detector response. To be able to 
generate signal events in practice, we split the processes 
into separable steps, in order to diminish the phase space 
from the otherwise up to ten-body final-state processes. 
In the first step the inert states are produced on shell, 
as in Eqs. ©-©. In the following steps, i.e. Eqs. (J9j)- 
(fTTJT> . the inert scalar particles are also taken to be on 




FIG. 4: Feynman diagrams contributing to pp — > 41 + fir in 
the IDM. 

shell while keeping the virtuality of the gauge bosons 
fully general. As a check of the validity of this approx- 
imation, we note that the inert particles' resonances for 
our benchmark models are indeed narrower than gauge 
bosons'. In all cases, the width of the ^4° is small, of the 
order 10~ 4 - 10" 5 GeV, due to the small mass difference 
to H . 8 The width of H± varies more, but is still smaller 
than the W width for all our models except one, which 
anyway has mass differences that allow both W and the 
inert state to be on shell simultaneously. Moreover, the 
most important contribution will come from direct pro- 
duction of A pairs, and the models for which the produc- 
tion of H gives a significant contribution to the signal 
coincides with large enough Am H ±^a = m H ± — m^a to 
allow A to be on shell in Eq. (|9|). 



A. Production of inert scalars via gauge fields 

In this subsection, we discuss the general expectations 
of the > 4 lepton signal strength from inert scalars pro- 
duced via gauge bosons. Some of the contributing di- 
agrams are shown in the first three panels of Figure [4] 
As the gauge couplings arc fixed, the production cross 
sections of the heavier inert states arc fully determined 
by their masses, and their decay patterns by their mass 
splittings: 

Amfj±H a = m H ± - (12) 
Am H ±Ao = m H ±—mAa, (13) 
Am^o H a = m A a—m H o. (14) 



This does not, however, make A sufficiently long lived to give rise 
to displaced vertices. 
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m A o[GeV] m A o[GeV] 



FIG. 5: Cross sections a into four or more leptons from 
H + H~ and H ± A° production via gauge bosons (in units of 
10 logarithms of a in fb). Here m H o = 70 GeV. This contri- 
bution to the cross section is independent of the Higgs mass. 



FIG. 6: Cross sections a into four or more leptons from A A 
and H + H~ production via the Higgs scalar h (in units of 10- 
logarithms of a in fb). In this plot m H o = 70 GeV, i4 = 
and m H ± — 220 GeV. Only Higgs production via gluon fusion 
is included. 



The processes in Eq. ©-([S]) can give rise to final 
states with four or more leptons if the mass hierarchy 
is m H a < m,A° < tuh ± , and this is the mass hierarchy 
we will consider in this section. For large Higgs masses 
( m h ^ 160 GeV), the EWPT constraints require F ± to 
be the heaviest state and the mentioned mass hierarchy 
is then just a consequence of H° being the DM candi- 
date. Apart from that, the Higgs mass has no impact on 
our results in this section. 

As we will see, even for optimal parameter values, the 
gauge-mediated contribution to a four-lepton signal in 
the IDM will not be enough to render the model de- 
tectable. Here we merely study under what conditions 
the contribution from gauge mediated production can be- 
come non-negligible, and we will turn to the more signif- 
icant contribution from gluon fusion in the next section. 

As our interest is in the detection of leptons, the 
branching ratios A° -> H°£ + £~ and iJ ± — » A /^ are 
important. For very small mass splittings Am^o^o the 
Bi(A° — > H°£ + £~) can be large, but give rise to leptons 
that are too soft to be isolated. For increased mass split- 
ting, decay modes into the more massive quarks open up, 
and the branching ratio into leptons decreases, approach- 
ing 6.7% which is the result for an on-shell Z boson. A 
small Aitiaoh also gives larger Br(i7 ± — > A°W ± ) as a 
large mass splitting would kinematically favor decay into 
H°; especially if W becomes on shell. Again a small mass 
shift becomes weighed against the ability to produce hard 
enough leptons for detection. For a fixed Am^o^o, in- 
creasing m,H± will typically increase Br(i7 ± — >A°), but 
at the cost of lowering the production cross section of 
heavier H . 



In Figure we show the cross section for the gauge 
mediated contribution to the production of four or more 
leptons. We calculate the tree-level cross sections with 
MadGraph/MadEvent and apply a corrective factor, a so- 
called K-factorjOf 1.2 to achieve agreement with the NLO 
results in [H,|90[. 



B. Production of inert scalars via SM Higgs 

The SM Higgs production at LHC is dominated by 
gluon fusion - dominantly induced by the loop of a top- 
quark coupled to the Higgs boson [9l|. The couplings 
of inert particles to the Higgs can then give a significant 
contribution to the production of four leptons through 
the processes 

gg -> h-> A°A°,H + H-. (15) 

In the A A channel one obtains four leptons in the 
final states independently of the values of mjj± and 
Br(iJ ± — > A ). This process is shown in the last dia- 
gram of Figure [4j The signal strength will, apart from 
m H o and m^o , also depend on and X H o . Unlike the 
processes considered in the previous section, the study of 
this process is strongly related to the SM Higgs search 
and to the search for DM in direct detection experiments. 
Given m#o and m^, direct detection data constrains the 
coupling X H a between H° and the Higgs boson, which 
for a given mass m^o also limits the size of the Higgs 
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Benchmark 


mh 


m H o 


m A o 


m H ± 


A 


\ mill 
^2 




A A o 


\ H ± 


0"W tot /3_body 


av -l~l/~lZ 


cr SI x o n « 






IDM-A1 


300 


72.0 


110 


210 


72 2 








0.22 


1.3 


1.5 / 1.5 


7.9 / 8.0 


4 ■ 10~ 5 


0.108 


IDM-A2 


500 


71.8 


110 


230 


72 2 


io- 7 


-10 -3 


0.22 


1.6 


1.5 / 1.5 


7.9 / 8.0 


5 • 10" 6 


0.111 


IDM-B1 


320 


77.5 


105 


152 


-1.6 • 10 4 


0.080 


0.70 


0.86 


1.3 


1.1/ 0.83 


2.8 / 2.5 


16 


0.110 


IDM-B2 


550 


76.0 


140 


220 


-6.0 • 10 4 


0.38 


2.1 


2.5 


3.4 


1.2 / 0.95 


3.7 / 3.7 


18 


0.113 


IDM-C1 


320 


91.0 


120 


190 


-6.5 ■ 10 4 


1.3 


2.3 


2.5 


3.2 


210 / 


15 / 62 


0.13 


0.00105 


IDM-C2 


280 


81.0 


130 


190 


-3.5 • 10 4 


0.50 


1.3 


1.7 


2.3 


16 / 14 


2.8 / 13 


8.1 


0.00979 


IDM-C3 


550 


92.0 


140 


230 


-1.6 • 10 5 


2.7 


5.4 


5.7 


7.2 


19/0 


0.38 / 3.6 


7.5 


0.0106 


IDM-C4 


550 


85.9 


140 


230 


-6.0 ■ 10 4 


0.38 


2.1 


2.5 


3.8 


23/0 


6.1 / 10 


1.4 


0.0110 



TABLE I: Benchmark models. Masses in units of GeV and A2 (not directly relevant) is taken to its, by vacuum stability, 
minimal allowed value A™ m . Annihilation cross sections, at relative impact velocity v — > 10 _3 c, are in units of 10 -26 cm 3 /s for 
o"Vtot,3-body and in units of 10 cm /s for <th 77i7 z. Spin-independent cross sections a in units of 10 cm . 



coupling to A 

2 2 

^A" — \ H a H t, . (lb) 

To generate gg — > A°A°,H + H~ events we make use 
of MadGraph/MadEvent's implementation of the Higgs ef- 
fective theory, where the Higgs boson couples directly to 
gluons. The effective coupling between the Higgs boson 
and the gluons depends on the Higgs mass and we match 
the cross sections obtained with MadGraph/MadEvent to 
the next-to-next-to-leading order results for Higgs pro- 
duction via gluon fusion in the SM [92| 9 . 

At the largest Higgs masses the vector boson fusion 
could also start to become relevant, but we are conserva- 
tive in the sense that we do not include such, or other sub- 
dominant, Higgs production contributions to our IDM 
signal. In Figure |6] we show the IDM cross section to 
four or more leptons by the Higgs mediated interactions 
in Eq. (HU). 

C. Background 

The requirement of leptons in the final state enables 
a signal to be extracted from the otherwise huge QCD 
background at hadron colliders. In order to simulate the 
SM background in the > 47 + fir channel, we include the 
following SM processes: 

VVV, ZZ, tiZ, ti, bbZ and titi, 

where V = W, Z are allowed to be off shell. 

Out of the contributions to VVV, WWZ is the most 
dominant contribution to our background and is the one 
we include in our analysis. We do not simulate WW 
processes, which are expected to be subdominant (93j . 



9 In fact, the default effective operator implementation in Mad- 
Graph/MadEvent-4.4.32 is not well suited for large Higgs masses. 
The deviation is as much as a factor 4.0, 4.6, 8.0 and 7.2 for 
mh = 300, 320, 500, 550, respectively, compared to the results in 
& 



We expect to be able to efficiently reduce these back- 
grounds in order to discriminate the signal: ZZ pro- 
duction is the dominant source of hadronically quiet Al 
events, but without invisible particles in the final states 
it can be efficiently removed by a cut on missing trans- 
verse energy. For IDMs producing leptons from off-shell 
Z bosons, the SM backgrounds including on-shell Z can 
be further discriminated by reconstructing the invariant 
mass of same-flavor, opposite-sign lepton pairs. The ttZ 
and titi backgrounds can also be reduced by vetoing b 
tagged jets, which should leave most of the IDM signal 
events. For the low background levels in the four-lepton 
channel, a significant contribution could come from fake 
leptons. This is difficult to properly take into account in 
a study based on Monte Carlo simulation, and should be 
estimated from experimental data. We comment further 
on this in our discussion of systematic uncertainties in 
Sec- IVIDl 

VI. ANALYSIS 

In order to study the signal expectations for IDM at 
the detector level, we define a set of benchmark models 
in Tabic HI The models are divided into three subsets: 

• The A models are constructed to test how strong 
the signal can be when inert states are produced via 
gauge interactions, and the direct detection signal 
will be very weak. These models do, however, give 
invisible Higgs branching ratios that are too low 
to pass the current LHC constraints on a heavy 
Higgs, and are therefore ruled out but kept here for 
illustration of the strength of the gauge-mediated 
production. 

• The B models represent IDM scenarios that explain 
all of the observed DM. They only pass all con- 
straints if we add the systematic uncertainties to 
the XENON-100 limits, as discussed in Sec. llVDl 

Models with a Higgs mass above 600 GeV (where 
no Higgs search limits have been presented) should 
more easily pass all constraints. Such models 
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Benchmark 








Br H ±_,. A o 


&4.1 


IDM-A1 


18.56 


54.00 


29.19 


0.191 


0.11 


IDM-A2 


13.36 


42.65 


22.68 


0.293 


0.12 


IDM-B1 


69.43 


123.8 


70.51 


0.071 


0.095 


IDM-B2 


16.23 


36.36 


19.12 


0.008 


0.003 


IDM-C1 


27.63 


62.44 


34.05 


0.013 


0.008 


IDM-C2 


27.20 


56.25 


30.47 


0.002 


0.001 


IDM-C3 


14.05 


32.71 


17.06 


0.079 


0.003 


IDM-C4 


13.77 


32.71 


17.05 


0.070 


0.022 



TABLE II: Cross sections for processes where the interaction is mediated via gauge bosons Z/'y or W , in units of fb. 



Benchmark 


a gg^A a A° 


a gg^H+H- 


&4l 


Br(ft -> 2H°+2A°+H+H-) 


IDM-A1 


88.25 


7.46 


0.40 


(0 + 0.84 + 0) = 0.84 % 


IDM-A2 


4.66 


138.0 


0.32 


(0 + 0.09 + 3.6) = 3.7% 


IDM-B1 


783.4 


1198 


4.1 


(5.9 + 7.7+13) = 27% 


IDM-B2 


194.6 


386.7 


0.90 


(4.3 + 5.7 + 15) = 25 % 


IDM-C1 


2844 


162 


13 


(32 + 30 + 0) = 62 % 


IDM-C2 


1981 


39.55 


9.0 


(27.4+ 19.4 + 0) = 47% 


IDM-C3 


483.0 


558.8 


2.5 


(15 + 16 + 28) =59% 


IDM-C4 


194.5 


366.1 


1.1 


(4.5 + 5.7+15) = 25% 



TABLE III: Cross sections for processes where the interaction is via the Higgs boson, in units of fb. 



should be able to give similar Al+ fir signal features 
as IDM-B2. However, a weaker signal is expected 
since the production cross section of inert states 
will be smaller once \ H o is adjusted to pass the di- 
rect detection constraint (unless we again allow for 
systematic uncertainties as for IDM-B2). 

The model IDM-B1 has a Higgs boson mass of 320 
GeV, motivated by the excess seen by the CMS 
experiment around this value. This possibility was, 
however, ruled out by the new ATLAS limits Q 
that were presented during the preparation of this 
manuscript (see Figure [5]). 



• The C models are illustrative examples of models 
that pass all constraints, but have a relic density 
that explains only a fraction of the observed total 
cold DM content. They are chosen such that IDM- 
Cl and IDM-C3 just pass the XENON-100 con- 
straint, but have some margin to the LHC Higgs 
bound. IDM-C2 and IDM-C4 instead just evade 
current Higgs searches at the LHC, but have larger 
margins to the XENON-100 limits. 



All the benchmark models pass all the other experimen- 
tal and theoretical constraints listed in Sec. IIIII 10 In our 
detector-level study, we take these as our representative 
IDM models for a tetralepton signature with a heavy SM- 
like Higgs boson. In Table HI1 and IIIII we list the models' 
properties relevant for the four-lepton signal. 

These models may well show up in upcoming data 
from XENON-100 and LHC. The expected performance 
of LHC is an integrated luminosity of up to ~ 15 fb -1 
collected by end of 2012 with an upgrade to 8 TeV for the 
rest of this year. The increase in sensitivity in the SM 
Higgs searches is about a factor 1.6 due to the integrated 
luminosity being 3 times lar ger and a factor about 1.2 
due to the increased energy [94j. A factor up to about 
\/2 could also come from combining the ATLAS and CMS 
data. This means that all our benchmark models, except 
possibly IDM-C3, should be reached by exclusion limits 
from LHC Higgs searches by the end of 2012. Detection 
of the Higgs bosons in any of our benchmark models, at 
the 5cr level, would however require more integrated lu- 
minosity, and such Higgs bosons would most likely not 
be revealed before the LHC run at 14 TeV. 

The cross section sensitivity of XENON-100 will also 
improve by an order of magnitude by the end of 2012 
[05L |9(| . This is enough to start to probe all our bench- 



The models IDM-C2, IDM-C3 and IDM-C4 give a 
relic DM contribution of 10% to f2cDM, and IDM- 
Cl gives 1% of ficDM- 



The high Higgs mass in combination with large couplings actually 
renders the IDM-C3 model marginally in violation of the, somewhat 
arbitrary, choice for the tree-level unitarity limit given in Sec. IIIII 
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mark models, except for the IDM-C1 and possibly IDM- 
C4 model (and of course the IDM-A models). The 
planned XENON-IT is expected to improve the sensi- 
tivity by more than an order of magnitude [95l l96j . 

We therefore consider the complementary four-lcpton 
plus missing energy channel as a potential step to further 
pin down or discover an IDM signal. 



A. Event generation 

We generate signal and background events with the 
MadGraph/MadEvent-4.4.32 package. From a user speci- 
fied process, MadGraph creates Feynman tree-level ampli- 
tudes (including effective operators and using a HELAS 
[97| implementation for the helicity amplitude calcula- 
tions) for all relevant hard subprocesses. Once events 
are generated with MadEvent they are passed to Pythia 
[8?| for hadronization and decay. The events are then 
passed to the Pretty Good Simulator PGS [H| to mimic 
the detector response. 

For each background and signal process, we generate 
events corresponding to an integrated luminosity of at 
least 10 times the integrated luminosity for which we 
make predictions. In a few cases, however, we were 
limited by computer power, and for the IDM-A mod- 
els and the SM backgrounds we have generated events 
corresponding to at least 3000 fb _1 , except for bbZ and 
tt production for which we have generated 220 fb _1 and 
160 fb _1 , respectively. 



B. Settings 

We consider proton-proton collisions at 14 TeV, us- 
ing the standard cteq6ll for the parton distribution func- 
tions (98|. In Pythia, we include initial- and final-state 
radiation but not multiple interactions. For our PGS 
settings we choose the options that mimic the ATLAS 
detector with a cluster finder cone size of AR = 0.4 
for jet reconstruction, and keep the other parameters as 
they are given by default in pgs_card_ATLAS.dat in Mad- 
Graph/MadEvent-4.4.32. 

For the cases where we generate events including jet 
matching (see S ec. IVID|) . we use the so-called MLM 
scheme |00l . llOOj with the minimum Kt jet measure for 
the phase space separation between partons set to 20 
GeV. 

The lepton isolation criteria are an important part of 
the lepton object definition in order to distinguish them 
from leptons that could have originated in jets. For elec- 
trons, PGS does this by default by requiring that the 
transverse calorimeter energy in a (3 x 3) cell grid around 
the electron, excluding the cell with the electron, has 
to be less than 10% of the electron's transverse energy 
and that the summed pr of tracks within a AR = 0.4 
cone around the electron, excluding the electron, is less 
than 5 GeV. To mimic the ATLAS detector response, 



we also ignore electrons with a pscudorapidity rj within 
1.37 < \r]\ < 1.52 jjxj. For muons, that are not isolated 
by default in PGS (and we do not make use of the clean- 
ing script that is the default in MadGraph/MadEvent), we 
require the summed px in a AR = 0.4 cone around them, 
excluding the muon itself, to be less than 10 GeV to de- 
fine them as isolated. For each lepton we also require a 
minimum distance of AR = 0.4 from the nearest lepton 
or jet (as reconstructed by PGS). 



C. Cuts 

In order discriminate an IDM signal from SM back- 
ground events, we perform cuts sequentially on the de- 
tector simulator's reconstructed particle data. 

To illustrate our cuts, we show in Figure [7] the event 
distributions after each cut. The plots include two of 
our benchmark models together with the total SM back- 
ground and two of its main subprocess contributions in 
this tetralcpton + fir channel. These are the cuts spe- 
cific to our IDM study: 

• First, we require four or more isolated leptons. In 
order to make lepton isolation and event triggering 
in the four-lepton channel robust, we will require a 
leading lepton with p l j, > 20 GeV and that each of 

the additional leptons have p%' ZA > 10 GeV. 

• In order to reduce the ZZ background efficiently, 
we require the missing transverse energy ($t) hi 
each event to be larger than 25 GeV, as illustrated 
in the (upper left) panel of Figure [7l 

• Wc reject events with any pair of same flavor and 
opposite sign (SF-OS) leptons among the > 4 lep- 
tons with an invariant mass that falls within the 
range of the Z resonance, 75 GeV< m[J < 
105 GeV. We refer to this as our Z veto. The (up- 
per right) panel in Figure [7] shows the distribution 
of events by the pair of SF-OS leptons giving an 
invariant mass closest to 91 GeV. 

• The tiZ background can be fairly efficiently dis- 
criminated against by requiring no b tagged jets in 
the event, as illustrated in the (bottom left) panel 
of Figure [7J Because of the displaced vertices from 
b quarks, the background from tt and bbZ could 
be further discriminated against by u sing a cut on 
the impact parameter for muons [10l| . Such an im- 
provement is beyond the scope of this paper, since 
it cannot be done within the standard PGS detector 
simulation that wc use. 



In the (bottom right) panel of Figure[?l we show the 
distribution of events in the minimal SF-OS dilep- 
ton invariant mass (minimal since each event has 
at least four leptons, and may contain more than 
one pair of SF-OS leptons). This invariant mass is 
expected to be low for our benchmark models, as 
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FIG. 7: Top left: Missing transverse energy distributions in events with four isolated leptons. Top right: The invariant mass 
of SF-OS lepton pairs after the cut on $t has been applied. Bottom left: Distribution of b tagged jets, after the cut on fir 
and Z veto. Bottom right: Invariant mass distribution for the SF-OS lepton pair producing the minimal such value per event, 
after all other cuts have been performed. The shaded grey regions indicate the cuts on each quantity. 



the Z decays off shell, and we require the minimal 
invariant mass to be < 50 GeV. 

For the signal events, the position of the peak in the SF- 
OS dilepton invariant mass distributions is slightly below 
the mass difference Am^o H o in a given model. The large 
fluctuations in the minimal invariant mass distribution of 
the total SM background (bottom right panel of Figurc[7]) 
come from the low statistics of our tt sample; only six tt 
events are left after the cuts, five of which lie in the 15-25 
GeV bins. This makes it difficult to say something about 
the distribution of this specific background contribution. 
What we can see is that if the tt events could be vetoed 
in some way, for example using the impact parameter 
for muons mentioned above, then the SF-OS dilepton 
invariant mass distribution can be used as a signature to 
clearly distinguish our models from the background. 



A characteristic of our benchmark models is that the 
signal leptons originate in off-shell Z bosons. Therefore 
our signal efficiency is sensitive to the isolation criteria 
and the minimum pr requirements on the leptons. Mod- 
els with larger Am^o^o, which allow A° to decay to on- 
shcll Z, would be more difficult to detect since in this case 
the signal cannot be distinguished from the background 
using the Z veto. 

D. Sources of systematic uncertainties 

Systematic uncertainties in our 41 + J£t signal study 
are due to limited statistics in some of our background 
samples, sensitivity to lepton efficiencies and fakc-lcpton 
contributions to the background. 

In our statistical analysis we fix the signal and back- 
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TABLE IV: Cross sections for backgrounds that require fake 
lcptons to produce four final-state leptons (the WWW,WZ 
processes). We also test if our backgrounds are sensitive to 
including jet matching. None of these effects seem important 
if compared to the values used in our final analysis presented 
in Table IVl (where jet matching is not included). We have 
required four isolated leptons and add the respective cuts for 
each column successively. Results are presented in units of 
10" 2 fb. We have included K factors of 1.6 for ZZ(j), 1.9 for 
the WZ{j) [ToJ and 1.4 for tiZ [j~03| . 



ground cross section expectations to our average results 
but, as mentioned in Sec. I VI Al generation of enough bbZ 
and tt events were limited by computer power. We trust 
that our cuts remove any contribution from bbZ, but the 
tt contribution gives an uncertainty in our background 
estimation. Our tt sample consists of only five events 
after all cuts, and for a Poisson distribution the upper 
expectation value is 9.3 events at 90% CL Taking this 
upper value as the average tt result instead would in- 
crease our total background cross section with only 30% 
(and a similar relative increase in the expected needed 
luminosities to discover the signals). 

The lepton efficiency is low for our models compared 
to the SM background. This is because the leptons in the 
model events originate in off-shell Z boson decay, and our 
signal predictions are thus sensitive to the lepton isola- 
tion and pt requirements. For comparison, if we use our 
plj, re quirements and decrease the lepton efficiencies, as in 
[l08j . both the signal and the SM background cross sec- 
tions are reduced by about 50%. Because of the increase 
of pile-up effects as the experiment reaches design lumi- 
nosity, the isolation criteria might have to be loosened 
and the px requirement raised in compensation. For our 
study, a raise to p^' 3A > 15 GeV would leave only 32% 
of the total background (completely remove the contri- 
bution from tt), while still leaving 20%-60% of the signal 
in our benchmark models. 

Since the background in this channel is low, we could 
be very sensitive to the contribution from fake leptons. 
In order to make use of PGS's ability to generate fake 
electrons, we show in Table ITVl the results of our cuts ap- 
plied on some SM processes that naively give three lep- 
ton final states (such as WZ) and include explicit jets. 
The table shows that these types of fake- lepton contribu- 
tions seem not to be very important. Likewise, we find 
that including jet matching would not alter the result in 
our final analysis that is presented in Table [V] (where jet 
matching, for consistency, is not included for neither the 
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TABLE V: Results of the SM background and total IDM sig- 
nal cross sections in units of 10~ 2 fb. We have required four 
isolated leptons and for each column, from left to right, we 
successively add our other cuts as desc ribe d in the text. K 
factors of 1 .6 for the ZZ background jl02f [ and 1.4 for ttZ 
and ti |103( ] have been applied, as well as the K factors for 
the signal processes as previously quoted. 

backgrounds nor the models). 

A proper inclusion of backgrounds involving fake lep- 
tons has to be b ased on experimental data. In a recent 
ATLAS analysis [l09j of the Al + $ T channel, the sys- 
tematic uncertainty due to differences in fake rate be- 
tween simulation and data was estimated to be around 
10% for the background processes tt and ttZ. They 
also find that the Z+jets give a significant contribu- 
tion to the background, potentially dominated by elec- 
tron Brcmsstrahlung in the detector material that subse- 
quently pair produce leptons. However, these events are 
found to contain ffi r of 2 0-60 GeV and hard jets, as can be 
seen in Figure 2 in |l09j ] . Requiring < 3 jets with px > 40 
GeV and optimizing the fir cut could potentially reject 
this background effectively without loss of more than ~10 
% of the signal events in our benchmark models. The un- 
certainties in the estimation of the Z+jets contribution 
are however large and an inclusion of this background is 
beyond the scope of our phenomenological study. 

Sources of systematic uncertainties will not be included 
in the following statistical analysis. 



E. Results 

In Table |V] we show the results after the signal and 
background events have been passed through the PGS 
detector simulation as we successively perform the cuts 
described in Sec. IVI CI 

To obtain a statistical measure for when our signal 
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Model IDM-A1 IDM-A2 IDM-B1 IDM-B2 IDM-C1 IDM-C2 IDM-C3 IDM-C4 

3a evidence, P obs = 50% (ftr 1 ) 810 300 64 64 19 3.8 20 50 

5a detection, P obs = 50% (fb^ 1 ) 2300 820 180 180 53 9.0 55 140 

95% CL exclusion, P obs = 50% (fb^ 1 ) 280 110 30 30 13 3.1 14 20 

TABLE VI: The expected integrated luminosities needed at 14 TeV for a 3a and 5a detection in the inert doublet benchmark 
models. Alternatively, the expected luminosity needed for a 95% CL exclusion of these benchmark models. 



could be observed or excluded, we assume the number 
of events to be Poisson distributed. The probability of 
observing or fewer events is then 

N f>n -B 

P(N-B) = J2 j— , (17) 

n=0 

given that the background expectation value B is the 
true mean. For a one-sided 3(5)cr detection, we take the 
probability (1— P(N; B)) of having this number of events 
or more due to a statistical background fluctuation to be 
less than 0.13% (2.9 x 10" 5 %). With a signal expecta- 
tion S, the probability to observe such an excess signal 
is 1 — P(N; S + B), and we request this probability P b s 
to be 50%. 

In Table I VII we show the prospects for when detection 
or exclusion of our benchmark models at the LHC will 
occur. The quoted integrated luminosities are for a 50% 
probability to have at least 3<r evidence, 5cr detection or 
95% CL exclusion of the models. 

Because of the sometimes low statistics needed to de- 
tect these models, the use of Poisson statistics should be 
more correct than e.g. the commonly used rule of thumb 
of a 5a discovery when S > max(5, 5yB). In Appendix 
[B]this and other commonly used statistical measures are 
compared. For the benchmark models with the strongest 
signal, and thus the lowest number of expected events at 
the time of a discovery, Poisson statistics lead to about a 
factor of two larger required integrated luminosity than a 
naive Gaussian approximation. In Appendix [Bj we also 
show that increasing the prospect from 50% to 90% prob- 
ability to find evidence for a signal can require up to a 
factor of two increase in the required integrated luminos- 
ity. 

VII. CONCLUSIONS 

We have investigated the status of the IDM in light of 
the results from DM direct detection searches by XENON 
and the Higgs searches at the LHC. These experimental 
results complement each other in constraining the viable 
parameter space. 

We first set out to study the IDM in the regime where 
the model both provides a DM candidate and includes a 
heavy Higgs boson, thereby possibly alleviating the LEP 
paradox. Considering the model's ability to evade the 
SM Higgs searches, we investigated the effect of imposing 
the bounds from direct detection assuming that the DM 



abundance is set by thermal freeze-out of the inert H° 
particle. 

In particular, the combination of constraints utilized in 
this work completely rules out the so-called 'new viable 
region' found in [56j where H° masses are in the range 
of 80-150 GeV. Moreover, we conclude that the ensemble 
of constraints are in conflict with the IDM for its whole 
viable cold DM mass range if the models shall also incor- 
porate the Higgs boson in the mass range 160-600 GeV. 
This conclusion can be avoided if either 1) the canonical 
experimental bounds can be relaxed or 2) the IDM does 
not account for all the DM. 

We investigate the prospects of detection/exclusion in 
the near future of models belonging to these types of 'es- 
cape' scenarios. Adding the systematic uncertainties to 
the observational constraints, and at the price of some 
fine-tuning, we found that we can still obtain IDMs that 
contain both a heavy Higgs boson (>500 GeV) and a 
good DM candidate. We also looked into the possibil- 
ity that IDM explains only a fraction of the universe's 
DM content, and thereby more easily evades current con- 
straints from both LHC and DM direct detection exper- 
iments. Some of these models can be efficiently probed 
by the foreseen data from XENON and LHC before the 
end of 2012. 

The potential detection of a heavy Higgs boson and/or 
a signal in direct DM detection experiments in the viable 
IDM DM mass range, although these would be striking 
features in favor of an IDM-like scenario, would not exclu- 
sively point to the IDM. A way to pin down the identity 
of the new physics further would be to compare different 
complementary channels. The prospects for detection of 
the IDM in the 14 TeV LHC data have been studied pre- 
viously for channels with two or three leptons, together 
with missing energy (29l . [301 ] . In this work, we have inves- 
tigated the possibility of a four-lepton plus missing en- 
ergy signature at the LHC coming from IDM. The models 
with a heavy Higgs boson that evade the current con- 
straints typically have large couplings between the inert 
states and the SM-like Higgs boson. As a result, the pro- 
duction of four-lepton final states via gluon fusion Higgs 
production becomes a particularly promising channel to 
track, and even discover, the IDM during the early runs 
at LHC's design center-of-mass collision energy. 

We find that in the four-lepton plus missing energy 
channel our benchmark points, where the inert particles 
are mainly produced via the Higgs boson, should show up 
early in the 14 TeV LHC run. Our models IDM-B1, IDM- 
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B2 and IDM-C1 to IDM-C4) should be seen at integrated 
luminosities of 3.8-64 fb" 1 (9-180 fb" 1 ) at the 3cr (5a) 
CL. We can note that the IDM benchmark points that 
were studied in the previous works [29|, H3] for the dilep- 
ton and trilepton channels, only one survives the cur- 
rent direct DM detection and SM Higgs searches. Nev- 
ertheless, according to these references, our benchmark 
points satisfy properties, such as favorable Am^ojyo, that 
should also render them detectable in the dilepton and 
trilepton channels at integrated luminosities of 100-300 
fb~ x . We thus conclude that, compiling recent experi- 
mental constraints, the IDM with a SM-like Higgs boson 
heavier than about 160 GeV could very well first show 
up in the tetralepton channel. 
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NOTE ADDED AFTER PUBLICATION 



Give n the lat est results from the LHC of a discovery of a resonance at ~ 126 GcV consistent with the SM Higgs 
boson [ll4 1 11 2| . we present in Tab. IVIll five IDM benchmark models representative of a SM-like Higgs mass around 
126 G eV (ll3l | . We have also checked that thes e are in agreement with the latest results of the XENON-100 experiment 
[l!5| . A more precise loop calculation in [ll6| typically renders 0-3/ > 10 _47 cm 2 , and should shift a few of our lowest 
asi to values close to that bound. The viable DM mass region is 50 < mjj ;$ 75 GeV with the updated constraints. 

The first two models give a non- negligible contribution to the Higgs width through h —> HqHq. None of the other 
models give a significant contribution to the invisible Higgs width. The third model is indeed chosen to have a 
negligible coupling to the Higgs, although the H° mass is slightly below half the mass of the Higgs, and the next 
three models have inert masses that are too large for the particles to be produced on the Higgs resonance. For the 
fourth model, the relic density is obtained through coannihilations resulting in a small a . For the fifth point, the 
annihilation cross section receives non-negligible contributions from both 2- and 3-body final states. The generation 
of the relic density for the last point is entirely driven through annihilations into WW, 2-body processes , although 
the zero velocity annihilation cross section is driven by 3-body processes. The first and the last benchmark points lie 
on the limit of the viable Hq mass range. 

Note that these types of models were not explicitly included in our four-leptons analysis. Inert scalars produced 
via a 126 GeV SM-like Higgs boson can only decay to A particles with masses below m^/2. Therefore, unless the 
Hq mass is below ~ 40 GeV, the final state leptons from Aq decays into Hq will be too soft to be detected as isolated 
leptons. However, an H° mass below at 40 GeV are too low to give the correct dark matter abundance and to be 
compatible with our other experimental constraints. Therefore these models never give a strong tetralepton signal. 
In the case in which inert particles are produced via the gauge bosons, the four lepton production is independent of 
the Higgs mass. The production is governed by electroweak couplings, and there is no resonance like in the case of 
production via the Higgs boson. We have already argued in Sec. V. A. that in such a scenario, the detection of a 
tetralepton signal at the LHC is unfeasible. Such processes could however eventually give rise to di- and tri lepton 
signals [29|, H(| . Our models are chosen in such a way as to be similar to the benchmark mod els in [2!| [3(| while 
still being consistent with all the most recent constraints included in our work. See also jll7j for three proposed 
benchmark models. 



m H o m A o m H ± /j 2 \ H a X A o X H ± c?W/3-bod y c« 77 / 7 z o~ SI Q H ah 2 Br(h — >inv) 

53.0 120 130 2100 0.023 0.39 0.47 0.097/0.0089 1.8/0.095 1.5 0.115 26% 

54.0 140 110 2500 0.013 0.54 0.31 0.056/0.016 2.3/0.17 0.50 0.107 11% 

60.0 160 160 3624 -7.6 • 10" 4 0.70 0.70 0.16/0.15 4.5/0.14 0.0015 0.110 0.02% 

65.0 72.9 120 4200 8.0 ■ 10" 4 0.036 0.32 0.40/0.38 5.7/3.1 0.0013 0.109 

65.0 120 150 3640 0.019 0.34 0.60 3.1/1.9 20/14 0.69 0.110 

75.5 130 98.0 6900 -0.038 0.32 0.086 1.0/0.91 4.5/3.8 2.1 0.104 

TABLE VII: Benchmark models with 771^=126.0 GeV. Masses are given in units of GeV. Annihilation cross sections, at 
relative impact velocity v — > 10~ 3 c, are in units of 10 -26 cm 3 /s for av tot , 3_bod y and in units of 10 -29 cm 3 /s for av lln z- 
Spin-independent cross sections ct si are in units of 10 -45 cm 2 . A2 = 0.01 for all models. 
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Appendix A: Naturalness 

The IDM serves as an explicit framework where a 
heavy Higgs boson, up to around 700 GcV, can be in- 
corporated and still be in agreement with EWPT. While 
this possibility is interesting in itself, it has also served 
as an additional motivation for the model. Indeed, a 
larger Higgs boson mass could alleviate the fine-tuning 
in the SM and make the model more natural by pushing 
the need for new divergence-canceling physics to higher 
energy scales (l7j . 

Raising the Higgs mass within the IDM does however 
not necessarily lead to improved naturalness as compared 
to the SM . The new inert scalars contribute with 
additional corrections to the SM-like Higgs mass, as well 
as exhibit quadratic divergences of their own. This can 
lead to increased overall fine-tuning although a larger 
Higgs mass naively renders it less sensitive to corrections 
from new physics at high energies. 

Let F 2 {pi) be a quantity that depends on some inde- 
pendent input parameters pi . The amount of fine-tuning 
in F 2 associated with pi can then be taken to be Ap. , 
defined by fTTol ] 



Cutoff scale = 1 .5 TeV 



5F 2 
"F 2 



= A F — 



Pi 



(Al) 



A model is said to be natural, up to an energy scale A, if 
the total amount of fine-tuning is sufficiently small. The 
exact upper limit on Ap 7 . in order for the quantity not to 
be considered to be fine-tuned is somewhat arbitrary. 

The scalar masses are the parameters that receive the 
dangerous quadratically ultraviolet-divergent contribu- 
tions. Using momentum cutoff regularization, the one- 
loop corrections to the scalar mass parameters fi 2 = 
fi 2 + SfM 2 can be written (as in [l9[ ) 



SM, max(A^) 
3 IDM, max(AV-i 




400 ln 500 
m [GeV] 



FIG. 8: Fine-tuning A, without RG effects, as a function 
of the SM Higgs mass. IDMs are represented by dark (red) 
marks for the benchmark models and light (green) ones for 
models in the scan. The SM, given a cutoff scale of 1.5 TeV, 
is represented by the (blue) solid line. The circles show the 
result for IDMs using Eq. IA4l and the crosses the result using 
Eq. IA5I The thick part of the solid (blue) line corresponds 
to the remaining Higgs mass window allowed within the SM. 
The dashed (blue) line is the SM result with RG running of 
the couplings included. 



momentum cutoffs of the loop contributions from fields 
associated with the SM doublet and the inert doublet, 
respectively. In our case the relevant fundamental pa- 
rameters are Af , A, & Pi- We will start by focusing on the 
Aj to assess the models sensitivity to physics at higher 
energy scales. 

Taking p, = Af for F 2 = (j, 2 , (j, 2 j Eq. IA1I implies 



64tt 



[-8A t 2 A t 2 + (3g 



and 



6fi 



64tt 2 



V 2 )A 2 



is 



■(2Aa 



[(3 ff 2 + 5 ' 2 )A 2 g + 8A 2 A 2 2 

+7^3 ^ 



- 8AiA 2 : 
A 4 )A 2 2 



A 4 )A 2 ! 



(A2) 



(A3) 



where At is the top Yukawa coupling, <?' and g are the 
U(l) and SU(2) gauge couplings and we have assumed 
independent cutoffs Aj. The loop contribution from in- 
ternal gauge fields arc sufficiently small that A lg will be 
irrelevant compared to A ( . For large scalar couplings 
the most relevant ones will be An 12 and A22 21 ~~ the 



dln^?2 



A' 11 ' 2 = - 

1 " din A? 



(A4) 



For each model, we take the fine-tuning to be 
A = max(| A^a 2 1)- Specifying an acceptable level of fine- 
tuning thus determines the cutoff scale up to which the 
theory is natural without introducing any new physics. 

In Figure [S] we plot the fine-tuning A for a given cutoff 
scale of Aj = 1.5 TeV. This cutoff scale corresponds to 
the perturbativity scale (in the SM) at which the one- 
loop RG corrections to the Higgs self-coupling grow to 
the same level as its tree-level value for a mu ~ 700 GeV. 
In [13] it was also used as the upper naturalness scale 11 , 



In 1171 the no-fine-tuning scale associated to the Higgs mass 
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Model: IDM-B1 




0.1 1 10 100 1000 



Q [TeV] 

FIG. 9: The running of IDM parameters with energy scale 
Q for one of our benchmark models (IDM-B1). Where the 
curves flatten out at around 150 TeV is when we terminate 
the calculation because the perturbativity limit of at least one 
\i > 47T is reached. The rr?„ Q and /j^ curves are normalized 
into units of their values m H o , Mi,o a t the scale Q = mu- 



and it was argued that with such a high scale one can 
no longer be certain that any new physics canceling the 
divergences will be observable at the LHC. 

The plot includes all the IDMs from our scans that 
give a relic density in accordance with WMAP and pass 
all the constraints in Sec. |TTT] except the Higgs bounds 
from LHC 12 . The solid blue line shows the results within 
the SM, and the mass range 115 to 129 GeV (the only 
span left for the SM Higgs given the current LHC limits) 
is marked as a thicker part of the solid blue line. This 
shows that A w 10 for the SM. The kink on the blue 
curve around 350 GeV is when the fine-tuning goes from 
being dominated by Aa 4 to being more sensitive to the 
Higgs cutoff Ah- We see that this measure A gives a 
large fraction of the IDMs (green circles) that are less 
fine-tuned than the SM (A « 10), but also many models 
that are not. 

A similar measure to Eq. lA4l was used in fl9| . but with 
the running of the parameters up to the cutoff scale also 
taken into account. With A = max(|A^2 2 |) and the RG 

equations deduced from [l04j - 4l06l | 13 . we find that with 
the fine-tuning condition A < 5 on /x 2 2 our benchmark 
models are natural up to cutoff scales A = 1.0 — 2.4 TeV. 
Figure [9] shows the running of the IDM parameters in the 
case of our benchmark model IDM-B1. For comparison, 
the SM is now left natural up to A = 1.2 TeV, with 
the SM Higgs mass is bound to be m h < 129 GeV [|. 



quadratic corrections was derived to be A = 1.3 TeV and to be 
independent of the Higgs mass. 

Here we also impose the constraint in Eq. 17 of reference [l7jl even 
though we note that this does not qualitatively change the result. 
All our renormalization conditions were set at Q = mj. 



Model: IDM-B1 




2 4 6 8 



A [TeV] 

FIG. 10: Fine-tuning measures, with RG effects included, for 
one of our benchmark models (IDM-B1). The notation for 
the plot legend is that Af = A* 2 , and \/J2~A% = A (as 
defined in the text). 

This measure leaves half of our benchmark models less 
fine-tuned than the SM. In Figure [8] we also added the 
SM result, for A = 1.5 TeV, when the RG running of 
couplings is included. 

In [H they used A = max (^(A^ 1 /) 2 + (A^ 2 ) 2 ) , 

where the contributions to A Pi associated with pi = Xi 
are also included. Here, the tuning with respect to A^ has 
no significant impact (on our benchmark models), but 
we comment further on this below. In Figure [10] we show 
this measure together with the individual contributions 
to the fine-tuning for our benchmark model IDM-B1. 

We note, however, that in the case of our the bench- 
mark models, the large quartic couplings are compen- 
sated by large negative values of /x 2 . to give small masses 
to the inert particles. This introduces an additional 
source of fine-tuning, even if each scale /if and /x 2 ,, as- 
sociated with the two Higgs doublets, is individually not 
severely tuned. This we can incorporate by introducing 
a fine-tuning measure on, e.g., the mass of the lightest 
inert particle 

A m H o = d\nm 2 H0 = 1 d(/4 - A g0i u 2 /2Ai) 

2 - din A 2 ~m% din A 2 [ ' 

(j% din fj% Affov 2 91n/i 2 
m 2 Ha <91nA 2 m 2 H0 <91nA 2 

This reflects better the increased fine-tuning in mod- 
els with high \ho and low m#o values. It also cures 
the artificial large fine-tuning that arrises when /x 2 goes 
through zero and drives A A2 H to infinity (any tuning 

around /i 2 . = is irrelevant as it then does not con- 
tribute to any of the inert particle masses). The resulting 
A = max(|A™ 2 HC |) are represented by crosses in Figure [5] 

(without RG improvement), as well as the solid black 
curve in Figure [TU] for IDM-B1 (including RG improve- 
ment). As our benchmark models come with rather large 
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Ai this measure typically leaves them less natural. Aa is 
less than 5 up to cutoff scales A = 0.4—1.4 TeV when in- 
cluding the RG evolution. With this fine-tuning measure, 
our benchmark models can thus hardly be considered to 
be less fine-tuned than the SM. 

We here also note that the sensitivity to variations in 
Pi = Xi could be significant already at tree-level. The 
tree-level contribution to 



\ m H0 _ 
A; — 



91ni 



<91n A, 



(A6) 



already gives A = max(|A, H ° |) ~ 6 — 25 for our bench- 
mark models, and is independent of A,. This type of fine- 
tuning is however not directly related to the unknown 
contributions beyond the cutoff scale, and would be ab- 
sent if we take our A; to be fixed and known parameters 
for each model. 



Note that Pd(N;X) are distribution functions, whereas 
Pl,2 are requested probabilities. 

It can be convenient to phrase the probabilities Pi 2 
in terms of a corresponding number m,2 of standard de- 
viations (n-a) for a one-sided normal distribution. We 
define such a correspondence by 



P 12 = i 
1,2 2 



1 + erf 



( Hi£ 



\V2 

where erf is the Gaussian error function 

2 



crf(x) 



dt e 



(B5) 



(B6) 



Equation IB5I thus defines what we refer to as an n-a ob- 
servation, independently of the type of distribution func- 
tion P D . For example 3(5)er corresponds to 1 — Pi, 2 = 
1.35 x 10~ 3 (2.87 x 10~ 7 ). 



Appendix B: Statistical measures 

It is desirable to have a statistical measure of the in- 
tegrated luminosity C expected to be needed to detect 
a signal with cross section as above a background with 
cross section as- 

We denote the probability of observing A" or fewer 
events from a distribution with expectation value X by 



P D (N;X), 



(Bl) 



where the index D distinguishes between different distri- 
butions. In the following, D = G and P to denote Gaus- 
sian and Poisson statistics, respectively. B = B(C) = 
asC and S = S(£) = a$C denote the expectation values 
of the number of background and signal events, respec- 
tively. 

To claim that an observation of A^bs events is an ex- 
cess, i.e. to reject the null hypothesis of a background 
expectation B, it has to lie outside the interval specified 
by the background model's Pi confidence level (CL). For 
a one-sided bound, this requires N ^ s > AT m i n (>C), where 
Nmin is the minimum integer number satisfying 



P D (AU n ;P) > P. 



(B2) 



For such a future observation to occur with a probability 
P2, when the underlying true scenario has an expecta- 
tion value S + B, it is required that N m i n also fulfills 
A^min < N max (C), where A^ max is the maximum number 
satisfying 



l-P„(N max ;S + B) >P 2 



(B3) 



For a given distribution function P D (N;X), the system 
of equations (|B2[) - (|B3|) can then be solved to find the 
smallest required integrated luminosity C that has an 
integer solution AT: 



AU n (£,Pi) < A^< N n 



c(AP 2 ) 



(B4) 



If the number A" of events is Poisson distributed, then 
the one-sided cumulative distribution function P D = P P 
can be expressed as 



P P (N;X) 



-y(N + i,x) 

T(N) ' 



(B7) 



where T and 7 are the ordinary and the lower incomplete 
gamma function respectively, 



for integ' 

r(AT) n>o 



N Yi 
i=0 



(B8) 



Strictly speaking, A" can only take integer values - as it 
represents the number of observed events - and in gen- 
eral one can therefore not replace the inequalities with 
equalities in Eqs. (|B2[) - (|B3[) and still find a solution. 
The analytical continuation (i.e. the gamma functions 
in Eq. (|B8jl ) can, however, be practical to have at hand, 
even though the final results should always derive from 
a solution with an integer N. 

If the number A^ of events is instead Gaussian dis- 
tributed, then P D = P G with 



P G (N;X) = ± 



1 + erf 



N-X 
V2a 



(B9) 



where we take a = \J X to coincide with a Poisson distri- 
bution for large X. In this case, Eq. (|B4[) can be written 
in a simple form. The Pi (n\-o) CL one-sided upper limit 
on the background being smaller than the P2 (ri2-a) CL 
one-sided lower limit on the signal plus background now 
reads 



B 



nVB <N< S + B- n 2 VS + B. (BIO) 



The expression for the required signal S can then be put 
into the following algebraic form (if we relax the require- 
ment of Af being an integer): 



S > ni 



n 2 



4P> + AVBnx + n\ 



.(Bll) 
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Model 


3a, g 


3a P 2 =90%, Q 


3a P 2 =90%, V 


5a, g 


5a P 2 =90%, g 


5a P 2 =90%, V 


Pi = P 2 = 95%, V 


IDM-A1 


720 


1600 


1700 


2000 


3400 


3600 


1000 


IDM-A2 


240 


590 


630 


680 


1200 


1300 


380 


IDM-B1 


44 


120 


140 


120 


240 


290 


88 


IDM-B2 


44 


120 


140 


120 


240 


290 


88 


IDM-C1 


10 a (19) 


36 


44 


28 


66 


90 


30 


IDM-C2 


0.95" (5.7) 


5.8 


8.3 


2.6 a 


9.3 


16 


5.2 


IDM-C3 


11° (19) 


38 


49 


30 


70 


96 


30 


IDM-C4 


33 


97 


110 


91 


190 


230 


70 



"Requiring a signal expectation of at least 5 events gives the value 
quoted in parentheses in the first column. 

TABLE VIII: Integrated luminosities C, in units of fb _1 , required to detect our benchmark models under different statistical 
measures. The 3(5)<r columns give the required integrated luminosity in order to observe a 3(5)cr evidence (discovery) with 
a probability P2 — 90%, when the number of event counts is assumed to be Gaussian, g, or Poisson, V , distributed. In the 
columns with no P2 value quoted, the commonly used criterion S > 3(5)\/B has been used (i.e the Gaussian approximation 
in Eq. IB12f) . In the last column we give the required luminosity to have a 95% probability to exclude the models with at least 
95% confidence. (See the text for further information.) 



Although N should be an integer, we will follow com- 
mon practice and leave out this additional requirement 
when we present results for Gaussian distributions in Ta- 
ble I VIH[ For m = n and n<i = 0, this gives the commonly 
used criterion for expecting an n-sigma detection 



S > tiVb, 



(B12) 



which corresponds to a probability P2 = 50% to observe 
the required iV b s from a Gaussian distribution that, in 
fact, also spans over negative A^bs- For ri\ = n 2 = n, 



Eq. (jBlip gives the sometimes seen criterion 107 



S > n 2 + 2nVB. 



(B13) 



From these equations the minimum C is easily derived 
by substituting B = ctb£ and S = asC. 

This defines our statistical measures to determine the 
expected integrated luminosity needed to observe a ri\- 
sigma detection with a probability P 2 . Equivalently, this 
formalism also gives the expected integrated luminosity 
needed to exclude the signal expectation S + B at the P2 
CL with a probability P\ . 



In Table IVIIII we present integrated luminosities re- 
quired to detect our benchmark models with different 
probabilities Pi. 2 under different assumed distribution 
functions Prj for the number of event counts. 
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